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PROGRAMMER'S GUIDE FOR SUBROUTINE PL0T3D 

Identification 

PL0T3D -A Subroutine Which Generates Three-Dimensional Hidden Line Displays 
Author -Larry Gales 

Date -May, 197S- Center for Quantitative Science in Forestry, Fisheries 

and Wildlife, University of Washington, Seattle, Washington 98195 

Purpose 

PL0T3D is a subroutine package which generates a variety of three-dimensional 
hidden line CalComp-tj'pcr displays. The package features multiple plots per page, 
with or without annotation, or annotation without plots; multiple pages; automatic 
scaling; flexible titling; two-dimensional interpolation over the image space; 
choice of size and location of plots on a page; choice of viewing angle and viewing 
distance; and various options which affect the appearance of the plot. Each page 
of plots fits on a standard S-1/2 x 11 inch page with margins of sufficient size to 
permit inclusion in three-ring binders. For a more detailed description of the 
purpose and usage of PL0T3D refer to its user's guide (Gales 1978). 

Usage 

PL0T3D communicates with the calling program through: (1) an argument list, 
(2) common blocks, (3) files, and (4) a set of file manipulation subroutines. 



4 



-2- 



• Argument List: 

PL0T3D is invoked by the following statement in the calling program: 
CALL QQPL3D(TLF, ERF, DTF, DTFl, NX, NY, 

LBOX, LHIDE, XMIN, XMAX, YMIN, YMAX, 

ZMIN, ZMAX, XRICH, YRICH, DFAULT, 

XVIEW, YVIEW, ZVIEW, XC, YC, 

XSIZE, YSIZE, PAGE, TITLE, LOCTIT, - 

PXC, PYC, ERR) J 
where QQPL3D is the main entry point in PL0T3D. All of the arguments, except 
PXC, PYC, and ERR are input arguments only and are unaffected by the operation 
of the program. ERR is an output argument which is set greater than zero if, 
and only if, PL0T3D detects an error. PXC and PYC are input/output arguments 
and contain the location of the lower left corner of the previous plot upon exit 
from PL0T3D. PXC and PYC should never be assigned values in the calling program 
except prior to the very first call, in which* case they should be set to zeroes. 
The types, restrictions, and descriptions of the arguments are as follows: 
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ARGUMENT LIST 



NAftE 


TYPE 


RESTRICT- 
IONS 


DESCRIPTION 


TLF 


Integer 


See file 
descriptions 


The unit number of a file written 
by the calling program which 
contains titles which annotate the 
plot. 


ERF 


Integer 


See file 
descriptions 


The unit number of a file written 
by PL0T3D which displays error 
messages. 


DTP 


Integer 


See file 
descriptions. 


The unit number of a binary file 
written by the calling program which 
contains the X, Y and Z coordinates 
of the image points to be plotted. 


DTFl 


Integer 


See file 
descriptions 


The unit number of a binary scratch 
file written by PL0T3D which contains 
the coordinates of points to be 
plotted. 


NX 
NY 


Integer 


2 < NX < 30 
2 _< NY < 30 


The number of grid cells in the 
image space along the x and y axes. 


LBOX 


Integer 


0, 1, 2, or 3 


The box control for subroutine PICTURE. 
This parameter controls the structure 
of the three-dimensional window, or 
box, which contains the plot. The 
options are: 

LBOX = 0: No box lines. 
LBOX = 1: Surface and box lines. 
LBOX = 2: Surface, box, and side 
bars. 

LBOX = 3: Surface, box, side bars 
and labels. 

Normally, LBOX = 3. 

■ 


LHIDE 


Integer 


0, 1, 2, -1, -2 


The hidden line control parameter for 
subroutine PICTURE. The options are: 






G 


LHIDE =0: No hidden lines r amoved. 
LHIDE = 1: Hidden lines of plot but 

not of box are removed. 
LHIDE = 2: All hidden lines are 

removed. 

LHIDE 0: Visible portions of the 

underside are also plotte 






Normally, LHIDE = 
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RESTRICT- 

NAME TYPE IONS DESCRIPTION 



XMIN 


Real 


> -10^^ 


XMIN, XMAX, YMIN, YMAX, ZMIN, and 


XMAX 






ZMAX define a three-dimensional 


YMIN 




^ O Q 

< 10^5 


wijidow which encloses all the data 


YMAX . 






in the image space to be displayed. 


ZMIN 






Data outside the window are not 


ZMAX 






shown. If one or more of the 








following conditions hold: 








XMIN >_ XMAX, YMIN >^ YMAX, ZMIN >^ ZMAX 








the computer considers the limits to 








be unspecified by the user and computes 








new limits for the axis or axes in 








question which are sufficient to 








include all the data along that axis or 








axes in the image space. 


XRICH 


Real 


> 0 


XRICH and YRICH are the Ax and Ay incre- 


YRICH 






ments used in two-dimensional inter- 








polation (enrichment). If either is zero. 








no enrichment takes place. The user should 








note that the values of XRICH and YRICH 








should be coordinated not only with the 








^ data hilt pI^o with t*hp sizp of the window 








set by XMIN, XMAX, YMIN, and YMAX. If 








either is too small, the computer will 








consume too much computer time, whereas 








large values of XRICH or YRICH will affect 








the appearance of the plot. Note that 








interpolation should only be applied to . ^ 








data which are correctly organized on the 








binary data file. 


DFAULT 


Real 


>^ -10^^ 


DFAULT is the default value assigned to 








all. cells in the image space. DFAULT 






< 10^3 


is usually set to zero. 


XVIEW 


Real 


1 XVIEW 1 


The point of view in space, in inches, from 


■YVIEW 




> XSIZE 


which the plot is seen by the user. 


ZVIEW 




1 YVIEW [ 


Normally, XVIEW and .YVIEW should be 






> YSIZE 


negative, indicating that the observer 






|ZVIEW| 


is behind the plot so that it appears in 






> XSIZE 


front of him, while ZVIEW should be positive. 






1 ZVIEW 1 


. indicating that the observer is above the 






> YSIZE 


plot looking down. The user may have to 








experiment with different values for these 








viewpoints in order to obtain the desired 




1 


view of the plot. 
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RESTRICT- 

NAME TYPE IONS ^DESCRIPTION 



XC 
YC 


Real 


0.5 < XC 
0,5 < YC 


< 
< 


8 

10,5 


The location on a page of the lower left 
corner of a box which encloses the plot. 
The distances are in inches from the lower 
left corner of the page. 


XSIZE 
YSIZE 


Real 


XSIZE > .1 
XC + XSIZE 
YSIZE > .1 
YC + YSIZE 


< 
< 


8 

10.5 


The size of a box, in inches, which is to 
contain the plot. 


PAGE 


Logical 


T or F 








If PAGE = T (true), then the 'current page is 
skipped and the current plot is written on the 
next page. If PAGE = F (false), the current 
plot is written on the current page. 


TITLE 


Logical 


T or F 








If TITLE = T (true), then the current plot 
consists of plot titles only, with no plot. 
If TITLE = F (false), then both plot and 
plot titles are written. Note that in the 
former case, the entire area defined by 
XSIZE, YSIZE is available for plot titles 
(or any other text), whereas in the latter 
case, the plot titles are either squeezed in 
near the top of the area defined by XSIZE, 
YSIZE, or just above it (see LOCTIT). 


LOCTIT 


Logical 


T or F 








If LOCTIT is T (true), then the plot titles 
will be contained wholly within the region 
defined by XC, YC, XSIZE, and YSIZE. If 
LOCTIT is F (false), the plot titles are 
placed above and outside this region. 


PXC 
PYC 


Integer 


None 
(set by 
PL0T3D) 








The location of the lower left corner of the 
previous plot (in inches). PXC, PYC are reset 
to zeroes whenever PAGE is true. 


ERR 


Integer 


None 
(set by 
PL0T3D) 








An output argument which is zero if PL0T3D 
detects no errors, or is greater than zero 
if it does detect errors. See the user*s 
guide for a detailed explanation of the error 
codes. 
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• Common blocks: 

PL0T3D uses blank common and four labeled common blocks named /QQPLl/, 
/QQPL2/, /QQPL3/, and /QQXYZ/. Blank common serves as temporary storage for the 
image space (F), a counter for the number of points mapped to each image cell 
(NP), and two temporary storage arrays (TEMPX, TEMPY) , and is structured as follows: 

COMMON// F(30,30), NP(30,30), TEMPX(30,30) , TEMPY(30,30) 

REAL F, NP, TEMPX, TEMPY 

Therefore, the calling program must reseirve 4 x 30 x 30 = 3600 (decimal) words of 
workspace at the start of blank common, e.g., 

COMMON// WSPACE(3600) 

REAL WSPACE 

Since PL0T3D uses this area for temporary storage only and does not preserve values 
between calls, the calling program may access and modify this region in blank common 
any way it chooses, except that values stored in the first 3600 words will be 
destroyed when PL0T3D is called. 

The three common blocks /QQPLl/, /QQPL2/, and /QQPL3/ are used only for 
internal operations in PL0T3D and can be ignored by the calling program. Common 
block /QQXYZ/, however, serves a vital role in facilitating binary input and 
output of data points for the image space both within the calling program and 
PL0T3D. /QQXYZ/ is structured as follows: 

COMMON /QQXYZ/ X, Y, Z, FINI, EOFT 

REAL X, Y, Z, FINI 

LOGICAL EOFT 

where X, Y, Z are the coordinates of one data point, FINI is the value assigned 

to the end-of-file indicator (-99999.0), and EOFT is a logical variable which is set 

true iff an end-of-file is read. Any binary data file passed to PL0T3D by a calling 
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program should be read, written, terminated, and rewound by the special file 
manipulation subroutines QQRXYZ, QQWXYZ, QQWEOF, and QQREW, respectively, contained 
in PL0T3D. The calling program must explicitly set FINI = -99999.0 before any of 
these routines is invoked; otherwise, FINI will be undefined and the results 
unpredic tab le • 

• Files: 

PL0T3D uses four files named TLF, ERF, DTF, and DTFl. TLF is the unit number 
of a formatted file which passes plot title information to PL0T3D. It must 
contain four card images, each of which is at least 66 characters long, which 
label the top of the plot. TLF is normally written by the calling program, although 
it may. reside on an external file. It is automatically rewound by PL0T3D at the 
start and end of execution. 

ERF is the unit number of a formatted file written by PL0T3D which displays 
any error messages. 

DTF is the unit number of a binary file which passes the x, y, z coordinates 
of plot image data points to PL0T3D. DTF is normally written by the calling 
program using subroutine QQWXYZ and must be terminated by an end-of-file written 
by QQWEOF. The order of points on DTF depends on whether enrichment (interpolation) 
is selected. If no enrichment is called for, the points may be ordered randomly - 
If two-dimensional enrichment is called for, then the points must form a sequence 
of triangles. For a more complete description of enrichment and its effects on 
ordering, refer to the user's guide. DTF is automatically rewound by PL0T3D at 
the start and end of execution. 
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DTFl is the unit number of a binary scratch file which is written 
and read by PL0T3D if enrichment is selected. DTFl is automatically rewound 
by PL0T3D at the start and end of execution. 





The characteristics of 


the files used by PL0T3D 


are summarized as follows: 


FILE 


AUTOMATIC 


READ BY 


WRITTEN BY 


UNIQUE UNIT 


NAME 


REWIND 


PL0T3D 


PLOT 3D 


NUMBER 


TLF 


Yes 


Yes 


No 


Yes 


ERF 


No 


No 


Yes 


Yes 


DTF 


Yes 


Yes 


No 


Yes 


DTFl 


Yes 


Yes 


Yes 


Yes 



The column labeled "AUTOMATIC REWIND" means that subroutine PL0T3D rewinds the 
file at the start of its execution and then rewinds it again just before it returns 
to the calling program. The column labeled "UNIQUE UNIT NUMBER" specifies whether 
or not different file names may reference the same unit number. Here, all files 
must be unique. 

Subroutine PL0T3D does not check the files for format errors nor does it 
check to see if file names reference valid unit numbers. These types of errors 
will generally trigger error messages and actions which are peculiar to a given 
computer installation. 

• File Manipulation Subroutines: 

PL0T3D contains four subroutines which manipulate the binary data files. These 
routines are invoked by PL0T3D and also by the calling program which prepares data 
for PL0T3D, and are always used in conjunction with the common block /QQXYZ/ which 
holds the x, y, z coordinates for a single point to be written or read. The 
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routines are: 



QQRXYZ (FILE ): Reads the x, y, z coordinates of one data point from 

FILE and stores the coordinates in /QQXYZ/, 
QQWXYZ (FILE ): Retrieves the x, y, z coordinates of a point stored 

in /QQXYZ/ and writes them on FILE, 
QQWEOF (FILE ): ' Writes the point (FINI, FINI, FINI) on FILE to indicate 

the end-of-file, 
QQREW (FILE ): Rewinds FILE. 

As an example of the use of these routines in the calling program, consider 
the following fragment of computer code which writes 1000 data points on unit 1, 
terminates unit 1 with an end-of-file, and then reads the points back in. Note 
FINI must be established before the routines are called and EOFT must be cleared 
before testing for an end-of-file. Also note that X= Y = Z = FINI = -99999.0 in 
common block /QQXYZ/ when the end-of-file is read. The code is as follows: 
... 

COMMON /QQXYZ/ X, Y, Z, FINI, EOFT 

C ^- SET-END-OF FILE INDICATOR 

FINI = -99999.0 

C 

C REWIND TAPE 1 

CALL QQREW (1 ) ■ , 

C 

C WRITE 1000 DATA POINTS ONTO TAPE 1 

DO 10 I = 1, 1000 
X - ... 
Y = ... 
Z = ••• 

CALL QQWXYZ (1 ) 
10 CONTINUE 

C 

C END FILE TAPE 1 

CALL QQWEOF (1 ) 

i2 
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C REWIND TAPE 1, AND PRESET EOFT TO FALSE 

CALL QQREWd ) 
EOFT = .FALSE. 

C 

C BACK THE 1000 POINTS 

DO 20 I = 1, 2000 
CALL QQRXyzd ) 
IF (EOFT) GOTO 30 

• • • ~ X 

• o • ™ Y 
... = Z 

20 CONTINUE 

C 

30 CONTINUE 
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Structure 

The overall structure of PL0T3D is displayed in the following flow chart: 



Yes 



Write error 
messages on 
ERF and set 
ERR > 0 
» 



•Set constants 

•Save the values of the arguments in 

common block /QQPLl/ 
•Initialize internal values 
•Set logic switches if windowing or 

enrichment is called for 
•Rewind TLF, read four title cards from 

it, and rewind TLF again 
•Check for errors in the argument list 



- \ Were errors found? 



No 



{ Is a new page called for? 



No 



3t 



Yes 



Move 


the 


plot to a new 


page 


and 


set PXC = 0, 


PYC = 


= 0 





Move to location XC,YC on the current page 
and set PXC = XC, PYC = YC 



Are titles only to be plotted? 



No 



Rewind DTF/ DTFl 



Is windowing called for? 



No 



Is interpolation (enriciiment) 
called for? 



No 



Read data from DTF or DTFl and plot 
it along with titles 



Yps 



Plot titles 



Yes 



Call QQWIND which 
determines the bound- 
aries of a window 
which encloses the 
data 



No [Were errors found? 
Yes 



Yes 



Call QQ2D which reads 
DTF and writes, inter- 
polated points onto 
DTFl 



No 



t/ere errors found? 



Yes 



I Rewind DTF> DTfF 



Return 



EKLC 



Subroutines 

The following is a list and brief description of all subroutines contained in 
PL0T3D in alphabetical order: 



QQ2D; 



QQCK: 



QQCSIZ: 



QQER: 



QQF: 



QQNCHR: 



QQPAK: 



Called when two-dimensional interpolation (env/ichment) is 
selected. It transfers all points from file DTF to DTFl and, 
after every third point, it calls QQTRI which adds large numbers 
of interpolated points to DTFl. QQ2D assumes that three 
consecutive points on DTF are the vertex points of a single 
triangle. 

Checks the input parameters NX,' NY, XSIZE, YSIZE, XC, and YC for 
errors in value. 

Calculates the size of characters and the distance between text 
lines, for the four title lines, so as to insure that they will 
properly fit in the space allocated within, or above, the plot. 
Writes out error messages on file ERF and sets ERR to an appropriate 
nonzero value. 

Reads data from a file and maps each data point within the x and y 
axes of the current window onto a cell in the image sp^ce F. It 
also accumulates the number of points mapped to each image cell, 
in the array NP, and tallies the number of z values which lie 
outside the ZMIN to ZMAX limits. 

Calculates the maximum number of characters in the title lines, 
minus all trailing blanks. 

Packr text consisting of one character per word into 10 characters 
per word (for the CDC 6400 computer). 
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QQPLOT: 

QQPLT2: 

QQPTTL: 

QQPTXT: 

QQREW: 
QQRXYZ: 



QQTRI: 



QQWEOF: 



QQWIND: 



QQI«YZ: 



Calls subroutine PICTURE which generates the actual CalConip plotter 
commands to produce a three-dimensional hidden line plot. 
Calls the CalComp subroutine PLOT to initialize the plotter so that 
the origin Is set at a specified point. 

Calls QQNCHR, QQCSIZ, and QQPAK to prepare the plot titles for 
plotting. 

Plots the titles prepared by QQPTTL by calling the CalComp plotter 
routine SYMBOL which generates the plotting instructions. 
Rewinds a file. 

Reads one x,y,z coordinate triple from the binary data file and stores 
it in common block /QQXYZ/. If X = Y = Z = FINI, it sets EOFT = 
.TRUE. 

Accepts the x, y, z coordinates of three vertex points of a 
triangular region and generates a series of equispaced linearly 
interpolated points along the plane defined by the three points. 
The points are written on DTFl. 

Writes the point (FINI, FINI, FINI) on a file. This point 
establishes the end-of-file. 

Computes the x, y, and z coordinates of a cubical window which 
encloses the data. I^ is called if: 

XMIN _> YMAX, or 

YMIN >^ YMAX, or 

ZMIN >^ ZMAX. 

If any one or more of the above conditions is encountered it 
computes new values for the x, y, or z axes which completely 
contain the data. 

Writes the x, y, z coordinates for the point currently stored in 
/QQXYZ/ on a file. ^ ^ 



External Subroutines 

The following is a list and brief description of all subroutines called by 
PL0T3D which are not contained in the source deck, 

PICTURE: This is a very large and complex plotting routine, adapted from 

the program developed by Melvin L, Prueitt of Los Alamos Scientific 
Laboratory, which generates commands for the CalComp plotter to 
produce three-dimensional hidden line displays. This routine is 
non-ANSI and rather machine dependent. 

PLOT; This is a CalComp plotter routine which is used to set the origin 

of a plot. 

PUTC: This is a small machine-dependent subroutine which packs one 

character-per-word text arrays into n characters-per-word, where n 
is the maximum number of characters which can be stored in a 
computer word (n equals 10 on the CDC 6400). 

SYMBOL; This is a CalComp plotter command which is used to plot characters. 

Coding Information * 

• Literals and constants: 

The literals used in PL0T3D can be divided into seven classes: 

1) the integers 0, 1, 2 and the real numbers 0.0 and 1.0 used as 
initial values or offsets; 

2) the integer 1 and the real number 1.0 used as an increment /decrement; 

3) the integers 0, 1, 2, and 3 and the real number 0.0 used in 
arithmetic tests; 

4) the integers 1, 2, and 3 used as subscripts; 
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5) the integer -3 and real number 0.0 used as arguments; 

6) ' the logical constants .TRUE, and .FALSE.; 

7) the integers 1 through 13 used as error numbers. 

All constants are assigned values in the "CONSTANTS" section of each routine 
and are described in the following table: 

CONSTANTS 



NAME VALUE SUB BLOCK DESCRIPTION 



NXD 


30 


QQPL3D 


/QQPL2/ 


Maximum dimension of 
along the x axis 


the 


image space 


NYD 


30 


QQPL3D 


/QQPL2/ 


Maximum dimension of 1 
along the y axis 


the 


image space 


PAGEW 


8,5 


QQPL3D 


/QQPL2/ 


Width of a plotter page, 


in inches 


EXTRME 


10^0 


QQPL3D 


/QQPL2/ 


An extremely large value 




EPS 


10--30 


QQPL3D 


/QQPL2/ 


An extremely small value 




MINXS 


0.1 


QQCK 




Minimum 


allowed 


value 


of 


XSIZE 


MAXXCS 


8.001 


QQCK 




Maximum 


allowed 


value 


of 


XC + XSIZE 


MINYS 


0.1 


QQCK 




Minimum 


allowed 


value 


of 


YSIZE 


MAXYCS 


10.501 


QQCK 




Maximum 


allowed 


value 


of 


YC + YSIZE 


MINXC 


0./f99 


QQCK 




Minimum 


allowed 


value 


of 


XC 


MAXXC 


8.001 


QQCK 




Maximum 


allowed 


value 


of 


XC 


MINYC 


0. 499 


QQCK 




Minimum 


allowed 


value 


of 


YC 


MAXYC 


10.501 


QQCK 




Maximum 


allowed 


value 


of 


YC 


STDCSZ 


0.27 


QQCSIZ 




Size of characters, in inches, in 
the titles for a standard size plot 


STDDYl 


0.36 


QQCSIZ 




Vertical distance, in inches, between 
title lines for a standard size plot 



18 



NAME 


VALUE 


SUB. 


BLOCK 


DESCRIPTION 


STDDY2 


1.4 


QQCSIZ 





Vertical distance^ in inches, between 
title lines on a standard size plot 
when titles only are plotted, without 
the plot itself 


STDSZ 


6.0 


QQCSIZ 


~~ ~ 


The size, in inches, of a standard 
size plot along the x and y axes 


ROUND 


1.5 


QQF 


___ 


A rounding factor used when locating 
the grid cell into which an x,- y, z 
point is stored 


BLANK 


IH 


QQNCHR 




The blank character 


LINE 


0 


QQPLOT 




A parameter used by subroutine PICTURE 


SC 


1.0 


QQPLOT 


6 


A parameter used by subroutine PICTURE 


STDSZ 


6.0 , 


1 

QQPLOT 




The size, in inches, of a standard 
size plot 


TITSZ 


4.2 


QQPLOT 




The size, in inches, of a standard 
size plot when titles only are plotted 


XCDSPL 


0.9 


QQPLOT 





The standard size of XC when titles 
onTv aTe to be Dlotted 


YCl 


0.0 


QQPLOT 




The YC displacement when titles only 
ai'e to be plotted 


NCL 


66 


QQPTTt 




The maximum number of characters in 
a title line' 


NLINE 


4 


QQPTTL 




The maximum number of title lines 



• Word Lengths: 

All values in PL0T3D are assumed to be stored in full length single precision 
words. x\ll alphanumeric values are stored one character per word, except text 
strings packed by QQPAK, which are stored 10 characters per word. 
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• Naming Conventions: 

All subroutines and common blocks start with the letters "QQ" In 
order to insure uniqueness. In addition, variables in a common block which 
are not used in a given subroutine are represented by dumn^y variables of the form 
ZZnn or ZZZn where n is a digit. These dummy variables may span more than one 
array or set of names. For example, assume subroutine SBl only makes use of 
the variable FINI in /QQXYZ/. Then /QQXYZ/ is declared in SBl as follows: 

FINI, ZZZ2 
FINI 

Limitations 

PL0T3D checks for 11 error conditions. If any of these conditions occurs, it 
outputs an appropriate error message and returns to the calling program. For a 
description of the error conditions and messages refer to the user^s guide. 

PL0T3D does not check to see if the files are correctly formatted or if unit 
numbers are valid. Such errors are left to the computer system. 

Extensions 

PL0T3D can be expensive to run if large numbers of data points are to be 
read, generated, or displayed. Much of the execution time is consumed in the input 
and output of data points using binary read/write operations (formatted read/write 
operations are even slower). If sufficient computer memory is available, however, 
one can drastically reduce this time by storing the data points directly in memory 

o 20 
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SUBROUTINE SBl 
• • • 

COMMON /QQXYZ/ ZZZ1(3), 
REAL ZZZl, 
LOGICAL ZZZ2 



and accessing them through simple memory reference operations. The four file 
manipulation routines discussed above were written with just this possibility in 
mind. It is only necessary to alter these routines in an appropriate manner in 
order to simulate the binary read/write operations with retrieve/store operations. 
The latter are logically identical to the former but are an order of magnitude 
faster. 

Computer Resources 

• Storage: 

The object deck for PL0T3D occupies 4142 (octal) words of storage when 
compiled under the Minnesota Fortran compiler on the CDC 6400 computer. To this 
must be added the 3600 (decimal) words in blank common used for the image space, 
plus buffer areas needed for the four files TLF, ERF, DTF, and DTFl. On the CDC 6400, 
this amounts to approximately 4242 + 7020 + 4000 = 17162 (octal) words of storage, 
some of which may be shared by the calling program.* 

• Execution Time: 

The execution time for PL0T3D depends primarily on the number of grid cells 
in the image space and the number of data points generated by two-dimensional 
interpolation. The following table gives the approximate execution times in CPU 
(central processing unit) seconds on the CDC 6400 computer as a function of grid 
size and number of data points: 




*Note that the CDC 6400 has up to four instructions per word. Hence, other 
computers may require substantially more space for the computer code. 
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PLOT 


GRID 


APPROX. NO. 




OPTION 


SIZE 


POINTS 


CPU SECONDS 


No interpolation 


30 by 30 


24 


4.0 


No interpolation 


5 by 5 


24 


0.42 ■ 


2D interpolation 


30 by 30 


1200 


6.72 


2D interpolation 


5 by 5 


1200 


3.53 



Sample Run 

The annotated listing on the next few pages illustrates the control, program, 
and input data cards for a sample run, along with the associated output. The 
output consists of a set of plots and annotation on one page which show the effects 
of various plot options applied to a single file of x, y, z coordinates. Note 
that the input data cards are processed by a free form input system (Gales and 
Anderson 1978; Anderson and Gales 1978). 
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XPL3D#Ct*5)0000#T35#P?. PLCT30 TEST 

ACCbUNT#XXXXXXXX,XXXXXX. 

COMMENT . 

COM ME NT. ♦♦♦♦*♦♦♦♦*♦ ♦*♦******♦**♦♦♦*♦**♦♦**♦**♦**♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦**'»♦ 
COMMENT.* THE APOVF CARDS TDEhTIFY THF JOB* SPECIFY THF MFMrjPY * 
COMMENT.* REOUIPEmENTS (50000 OCTAL)» ESTIMATE THF CENTRAL PPHCFSSCR * 
COMMENT.* TIME (35 SE C ONOS ) » . A N D SPECIFY THE USERS ACCOUNT * 
CGMME NT. ♦**********♦**♦*♦***♦*******♦*♦*♦**♦♦♦*****♦ *♦**♦♦****♦*♦*♦*♦♦♦♦ 
COMMENT.' 
MNF»E«C#R»7# J. 
CCM^'ENT. 

COMMENT .**♦*♦♦*****♦**♦**♦****♦♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦***♦)*:♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦* 

COMMENT.* THF ABCV'E CARD CALLS THE MINNESOTA FORTRAN CrMPTL^P TC * 
COMMENT.* COMPILE THF EXECUTION PROGRAM FOR PL0T3D ♦ 
CCMMENT. ***************************************************** *********** 
CCMMENT . 

ATTACH, BP L 30, I D» BP L 3D . 
COMMENT. 

COMMENT. *********************♦****♦*♦**♦♦♦**♦♦♦♦♦*♦♦♦*♦♦♦**♦♦♦ 
COMMENT.* THE APOVF CARD ATTACHES THE PLOTTING ROUTINE TN BINARY FCFM * 
COMMENT.**************************************************************** 
COMMENT . 

ATTACH, BEE, ID-BFF. 
COMMENT. 

COMMENT. ***************************************************** *********** 
COMMENT,* THE ABOV^ CARD ATTACHES TKE FRFE-FOPMAT INPUT ROUTINES IN * 
COMMENT.* PINAPY FORM. * 
COMMENT ,*********************** **********^**********^******************* 
COMMENT. 
ATTACH, GRAFTN. 
L ieRARY,GRAFTN. 
COMMENT. 

COMMENT. *****'«'******************** ************** 

COMMENT.* THE ABOVE CARDS ATTACH THE LIPRAPY * 

COMMENT.* ROUTINES FOR PL3D AND DESIGNATE THAT * 

COMMENT.* FILE AS THE CURRENT LIBRARY. * 

COMMENT.**************************************** 

COMMENT. 

LOAD,BFF . 

L0AD,BPL30. 

LOAD, LGO. 

EXECUTE, XPL3D. 

COMMENT. 

COMMENT ************** *************************************************** 
COMMENT.* THF ABOVE CARDS LOAD THF EXECUTION, INPUT AWP PLOT ROUTINFS * 
COMMENT.* AND CAUSE CONTROL TO BE PASSED TO XPL3D FOR FXECUTION * 
COMMENT. ********************************************************** ******* 
CATALOG, TAPE99, RP»?. 
COMMENT. 

COMMENT. ************* ******************* ******** 
COMMENT.* THE ABOVE CARD SAVES THE PLCTTINf * 
COMMENT.* TAPE»-TAPE99-, SO THAT THF PLOTS MAY * 
COMMENT.* BE PREVIEWED BEFORE BPING SENT TO * 
COMMENT.* THE PLOT PUEUE. * 
COMMENT.**************************************** 
COMMENT. ' ' 

COMMENT.************* ******♦*♦♦**♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦.♦,.« 

COMMENT.* THE BELOW CARDS ARE THE EXECUTION PROGRAM XPl3n » 
COMMENT. ********************************************** ************ ****** 
COMMENT. 
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♦ EOR 



POCGPAM XPl3n (T KPUT.nUTPUT» TAPFl. TAPF2. TAPF3* TAPF'^-TMPUT, 
, TAPFfr-OLTPUT) 



PURPOSE 

DFMPNSTPATES TH F USE OF THF PLCTTINC- PACKAGF PL'^T3D 



C 
C 
C 
C 

c 

r 

C-IDFNTIFICATinN 
C 
C 

c 
c 
c 

c 

r 

c 
c 
c 

C-GLOBAL VARIABLES 
r 

C 

c 
c 

r 
C 

c 

c 



DATE « 
LGCATinN i 



LANGUAGE t 
INSTALLATICNJ 
PRTCPAhiMFRS : 



AUGUST 1977 

CENTER FCP QUANTITATIVE SCIENCE 

TN FORESTRIES* FISHERIES, AND V'TLOLTFF 
UNIVERSITY OF WASHINGTON 
SEATTLE, WASHINGTON 98195 
ANSI STANDARD FORTRAN 

CDC 6400/CYRFP 73 UNDER SCOPE 3. A H.S. 
L. GALES, BAILEY, J. SANISLO 



ERIC 



c 
c 

c 

c 
c 

C 

c 
c 
c 
c 
c 
c 
c 
c 

C 

c 



************************** ********'*** 

* BLANK COMMON ♦ 

* USED FOP STORAGF OF IN^UT VARIAPL^^S RFAD * 

* BY INPUT PACKAGE. * 

* -WSPACE- IS USED 8Y THE PLOT ROUTTNFS. * 

:i,***^**************************************-^* 



COMMON / / 


WSPACE(36CC), 










• 


N y. 


NY, 


LBOX, 


LHIDE, 


XMTN', 


XMAX « 


• 


Y MIN, 


Yr^AX, 


ZMIN, 


ZMAX, 


XPTCM, 


YOICH, 


• 


DFAULT, 


XVIFW, 


YVIEW, 


ZVIEW, 


xr, 


YC» 


• 


X SIZE, 


YSIZE. 


PAGE, 


TTTLF, 


LnCTIT 




REAL 


WSPACE 












INTEGER 


NX* 


NY» 


LBTX, 


LHIOE 






REAL 


X MIN, 


XMAXf 


YMIN, 


YMAX, 


ZMIN, 


ZMAX, 


• 


XRICH, 


YPICH, 


DFAULT, 


XVIFW» 


YVTFW. 


7VTFW, 




XC, 


YC 










LOGICAL 


PAGE, 


TITLE, 


LCCTIT 









OFF INITIRNS 

WSPACt(*) • 
NX, NY = 
LBOX « 
LHIDt • 
XMIK, XMAX, 
YMIN, YMAX, 
ZMIN, ZMAX » 

XRICH, YRICH 

OFAULT » 
XVIEW,Y\/If W» 
ZVTEW = 

XC.YC » 



STORAGE AREA FOR PLOT PQUTINtS 

NIJM3ER GF X AND Y TFLLS IN THE IMARF "^^Arc 

ROX CONTROL FOR SUiJROUTINE -PICTUPF- 

HIDOEN LINE RFMOVAL FLAG FOR SUBPniiTINF -pyCTLRt- 



UEFINE A RECTANGULAR P AP ' L L E L 0 P I P F D .WH TC H cmclCSF'^ 

TriF DATA TO BF PLOT.cD 
THE DELTA X AND DELTA Y INCREMENTS HSFD IN TWf- 

DIMPNSIONAL ENRICHMENT 
THE DEFAULT VALUES ASSIGNED TO ALL POINTS 

THE POINT CF VIEW FROM WHICH THE PLOT TS n^^SFPVFP 
(IN INCHES) 

LOCATION ON PAGE OF LOWER L'^FT CGRNFP OF PLOT 
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C 
C 
C 
C 
C 
C 

c 
c 
c 
c 
■c 
c 
c 



c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 



c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 



c 
c 
c 
c 
c 
c 

c- 

c 



(IN INCHES) 

XSIZE#YSIZE » PHYSICAL ^IZP OF PLOT (IN INCHES) 

PAGE « .TRUt. IFF THF CURRENT PLOT IS TCI /^PPfAR ON th(- N'f^XT 

PAGF 

TITLE - .TRUE. IFF THf CURRENT PLOT rnNSI*^T«^ OF ONLY TFXT 

* * 

* DATA COMMON PL3D1 * 

* USED FOR STORAGE OF THE DATA POINTS * 

* ♦ 

C0MM0N/PL3D1/ X(6)» Y(8)» Z(fl) 
REAL X# Y# Z 

DEFINITIONS 

X(*),Y(*)*Z(*) - THE LOCATION OF EACH DATA POINT 



* * 

* VERTEX COMMON PL3D2 ♦ 

* CONTAINS THE LOCATION OF THE VERTICES * 

* FOR EACH TRIANGULAR REGION ♦ 

* ♦ 



C0MMCN/PL3D2/ EX(3,3)»£Y(fi»3)»EZ(P»3) 
REAL EX> EY> EZ 



DEF INITIGNS 

EX(*,*) • 
EY(*>*) - 
EZ(*>*) - 



THE X COORDINATE OF EACH VERTEX 
THE Y COORDINATE OF EACH VERTEX 
THE Z COORDINATE OF EACH VERTEX 



* * 

* COMMON QQXYZ ♦ 

* USED IN QQPL3D FOR ACCESSING THF DATft * 

* POINTS AND CHECKING FOR AN END-OF-FTLE * 

* * 

**4i4: ********************************** ******* 



COMMON/QQXYZ/ Z Z Z 1 ( 3 ) » F I N I , 
ZZZ2 

REAL ZZZ1> FINI 

LOGICAL ' ZZZ2 



DEFINITIONS 
FIM ■ . 



-cqqsg.O, AN ENO-Op-FILE MARK 



ERIC 



•LOCAL VARIABLES 

INTEGER TLF, 



EPF 



DTFj 



DTFl » 



INF, 



ECF# 
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C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 

r 
C 
C 
C 
C 
C 
C 
C 
C 
C 

c 
c 
c 
c 
c 
c 

r 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 



RtAL 

DEFINITIONS 

TLF • 

fPF » 
DTF » 
DTFl » 
INF » 
ECF « 
DCF « 
ERP = 
C(*) « 
I, J,K » 
PXC>PYC « 



OCF, 
PXC, 



EKP, 
PYC 



C(4?>), I» 



J, 



-SUBROUTINES- 



-FXTERNAL^- 



QCPL3D 
PLOT, PLOTS 



OOWEOF 
QORFAD 
•FILES 



TLF 

FRF 

DTF 

OTFl 

INF 

DCF 

ecF 



-CONSTANTS" 



C 
C 

C-INITIALIZATIDNS- 

C 

C 



TITLE FILfc CONTAINING FOUR TITLES WHICM LABi^L PLTT 

THF ERROR MESSAGE FILE 

DATA FILE OF X, Y> AND Z COORDINATFS 

SCRATCH DATA FILE ' 

TNPUT FILF CONTAINING VARIABLE VALUES 
FILE USED TO ECHO INPUT 

DECLARATION FILE USED TO DECLARE INPUT VAPIAOLES 
ERROR NUMBER 
THE TITLE TEXT 
LOOP INDICES 

LOCATION CF THE LOWER LEFT CDPNER OF TH? PPFVITLS 
"LOT (IN INCHES) 



EL 

NMLIST 



THE PLOTTING ROUTINE 

THE INITIALIZATION AND TERMINATION PQUTTNES FTP 
THE PLOTTER (LOCATED IN THE GRAPHIC*; LIBRARY) 
WRITES AN END-OF-FILE MARK ON THE DFSIGNATFD FILE 
(LOCATED IN THE PLOTTING PACKAGE) 
READS INPl'T VARIARLFS IN FREE FORMAT 



INPUT, 

OUTPUT, 

INPUT* 

I N/OUT, 

I NPUT, 

OUT/IN, 

OUTPUT, 



FCRMATTED, 

FORMATTED, 

BINARY, 

BINARY, 

FORMATTED, 

FORMATTED, 

FORMATTED, 



textfile 

error messagp frp 
data file of c dor 0 in a t f 
scratch tata frf 
input file 
oeclap'ation file 
echo file 



DATA 


X / 


1.0, 3.0, 


2 .0, 


3.0, A. 2. 'j.n. 


4.5. 3 . C / 


DATA 


Y / 




2.0, 


2.4, 2.e» • 3.0, 


?.5, l.C / 


DATA 


TLF, 


ERF, 


DTF, 


DTFl, INF, 


CCF / 


• 


3, 


6, 


1, 


2, 5, 


6 / 


• 


DCF, 


FINI 


/ 






• 


' 1, 


-9<i999.0 


1 






• 


C / 


8*1HX, 


8+lHY, 


32*]HT / 





ERIC 



THE 

CALL PLOT? 



CALL TO -PLOTS- TNITIALI7F'; THE PLCTTFR 

2« 
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C 
C 
C 



PXC » 0.0 
PYC « 0.0 

/* THt CALL TO -NMIIST- DFCLA'?F'^ TNPl'T VARTAnLtS 
CALL NMLIST (3600# IH #nCF»ERF, ER" ) 

/♦ THE CALL TP -^COPEAD- RFADS INITIAL VALl'^S 
CALL QQREAD (INF,ECF»FPF» ERR) 



C-START- 
C 



DO 70i 

GaT0(I0#?0#30#'^0)#I 



10 
702 



DO 702 j-i»e 

Z(J) - 1.0 

CONTINUE 

/* READ PLOT PARAMETERS 
CALL OQREAD(INF,ECF»ERF, ERR) 
GOTO 999 
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Z(2) « 9.0 
Z(A) • 9.C 
Z(5) - 9.0 

/♦ READ PLOT PARAMETERS 
CALL OQREAD(INF,£CF, tRF, ERR) 
GOTG 999 



3C 
703 



C 

C 



AO 



DO 703 J«l>8 

Z(J) - 9.0 

CONTINUE 
Z(5) - C.C 

/♦ READ PLOT PARAMETERS 
CALL OOREAD(INF#FCF#ERF, ERR) 
GOTO 999 

/* nEAD PLOT PARAMETERS 
CALL OQREAD(INF,ECF,ERP» ERR) 
GOTO 998 



999 



705 
70A 

998 

1 FORMAT 



ERIC 



CALL FL(1> 1>3»2) 
CALL eL(2> ?#3»<i) 
CALL EL(3> 2»A,^) 
CALL EL(A, Z,5tb) 
CALL FL(5# 5»6#7) 
CALL EL (6, 5, 8,7) 
CALL EL(7# dt^itB) 
CALL EL(B» A,3»e) 
DO 70^ J-l#8 

DC 705 K»l>3 
WPITEOTF) 

CrNTINUE 
CONTINUE 

CALL CQWEDF(DTF ) 
WRITE(TLF#1) C 
( 8AI0 ) 

CALL C0PL3D(TLF» ERF,DTF, DTF1*NY,MY. L8 0X,LHTr.F# 
XMTN,XMAX,YMIN.YMAy.ZMIN, 7MAX, 
XRICH,YRICH» DFAULT,XVTFW» YVIEW» 
07 ZVIEWfXC* YC*XSIZE# YS HE >PAGf# TITLE* 

* LOCTIT* PXC*PYC» FRP ) 



EX(J,K),EY( J.K),FZ(J»K ) 
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c 
c 
c 



c 

r 
C- 

c 
c 
c 
c 

C' 

c 

c 
c 
c 
c 

c 
c 

C" 

c 
c 
c 



c 
c 



701 CONTINUE 

/♦ THE CALL TO -PLOT ( 0# 0# Q99 )- TFRMIN'ATFS Tt-f 
/* PLCTTFP 
CALL PL0T(0#0#999) 

STOP 
END 

SUBROUTINE FL(I* I1#I2>I3) 



•PUPPOSt 

CALCULATE? THE VERTICES OF EACH TRIANGULAR PEGIPN 

■ARGUMENTS 

INTEGER I» .Il> I2> 13 

DEFINITIONS 

I»I1»I2,I3 - INDICES" FOR THF CALULATIGN OF THE VERTKES 

•GLOBAL VARIABLE*^ 

note: FOR DEFINITIONS REFER TO Mf\IN PROGRAM 



REAL 



REAL 



C 
C 

C-START- 
C 



X (P). 


Y(B ), 


Id) 


X , 


Y, 


Z 


EX(8>3), 


EY(8, 3), 


EZ(P»3) 


EX, 


EY, 


EZ 


EX{ 1,1) 


» X(I1) 




EX(I,2) 


» X(I2) 




EX( Ii>3) 


» X{13) 




E Y{ I, 1) 


» Ydl) 




EY{I,2) 


» Y(I2) 




EY(I,3) 


« Y(I3) 




EZ( 1,1) 


» Z(I1) 




FZ( 1,2) 


» Z(T?) 




EZ( 1,3) 


» Z(I3) 




R ETURN 






FNO 







c 

C 

C- 

C 

c 
c 

C 

c- 
c 



SUeftPUTINE NMLI STdNDX, NAME, DCF, ERF, ERR) 

PURPOSF ; 

WRITE DECLARATION FILE USED 3Y FRFF-FORMAT INPUT ROUTINES 



•ARGUMENTS 

INTfeGFR 



INDX, 



NA^'S, 



OCF, 

•^8 



ERF, 



eoR 
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DEFINITIONS 

INDX - 
N&HF « 
DCF • 
EPF « 
ERR - 



C 
C 

c 
c 
c 
c 
c 
c 
c 
c 

C-EXTEPN6LS- 
C 

C OOINTL 

C 

C 

C-START 



1 



FGPMAK 
.e2H 
./62H 
./62H 
. /62H 
./62H 
./62H 
./62H 
./62H 



STAPT OF Variable stopage in blank comhtn 

I /□ LIST IDENTIFIER 
DECLARATION FILE 
ERROR FILE 

ERROR FLAG {REMAINS 0 IF NO ERROR*^) 



FREE-FORMAT INPUT ROUTINE TO PPOCP'^'^ VARJABLFS LIST 



REWIND DCF 
WRITE(DCF>1) 
REWIND DCF 

CALL QQINTL (INDX»NAMP,DCF,ERF> 
REWIND DCF 
i? fcTURN 



INTEGER 
REAL 



LOGIC AL 

S 

S 



NX, NY. LBDX, 

XMIM, XMAX, YMIN, 

Zf'Ax, XRICH, YRICH, 

WlEW* ZVTEW, XC» 
YSIZE 

PAGE, TITLE, LOCTIT 



EPP ) 



LMinE 

YMAy» 7MTN, 

DFAULT» XVTFW, 

vr, X5IZF, 



♦ EOR 



END 



On 
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2.0, YSIZE • Z.Ot 
$ 



/♦ 

/♦ INITIAL PLOT PARAMETERS */ 
/♦ 

■ XMIN • 1«0# XMAX » 5.0# YMIN • 1.0# Yt^AX ■ C,0, 
ZMIN » 0.0> ZMAX • 15. 0# XRICH ■ 0.05» YRICH • O.O''* O^AULT » 
XVIEW • 5,00 YVTEW » I3.22e# NX « 20# NY » 20» 
ZVIEW • 30.0, 

LBOX » 3, LHIPE » 2 $ 
/♦ 

/* pAPAf^ETERS FCP PLOT 1 */ 
/♦ 

XC • 1.5, YC » 6.0, XSIZF 
TITLE • .F., LOCTIT • .F. 

/* 

/* PAPAfETEPS FOP PLOT 2 */ 
/* 

XC • 5.0, LOCTIT « .T. $ 

/* 

/♦ PARAMtTEPS FOR PLOT 3 */ 
/♦ 

XC " 3.0, YC - 3.5, XSIZE 

/♦ 

/♦ PARAflETEPS FCP PLOT * /' 
/* 

XC • 3.0, YC • 1.*^, XSIZF » 5.0, 
TITLE = .T. $ 

*EOR 
*EOF 



CO, 



» 3.0, YSIZE » 3.0» LnCTIT ■ .F. 



YSIZr » 1.0, 
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